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ABSTRACT 



Context. The structure of the solar chromosphere is believed to be governed by magnetic fields, even in quiet-Sun regions that have 
a relatively weak photospheric field. During the past decade inversion methods have emerged as powerful tools for analyzing the 
chromosphere of active regions. The applicability of inversions to infer the stratification of the physical conditions in a dynamic 3D 
solar chromosphere has not yet been studied in detail. 

Aims. This study aims to establish the diagnostic capabilities of non-local thermodynamical equilibrium (NLTE) inversion techniques 
of Stokes profiles induced by the Zeeman effect in the Can /18542 A line. 

Methods. We computed the Can atomic level populations in a snapshot from a 3D radiation-MHD simulation of the quiet solar 
atmosphere in non-LTE using the 3D radiative transfer code Multi3d. These populations were used to compute synthetic full-Stokes 
profiles in the Can /18542 A line using 1.5D radiative transfer and the inversion code Nicole. The profiles were then spectrally de- 
graded to account for finite filter width and Gaussian noise was added to account for finite photon flux. These profiles were inverted 
using Nicole and the results were compared with the original model atmosphere. 

Results. Our NLTE inversions applied to quiet-Sun synthetic observations provide reasonably good estimates of the chromospheric 
magnetic field, line-of-sight velocities and somewhat less accurate, but still very useful, estimates of the temperature. Three dimen- 
sional scattering of photons cause cool pockets in the chromosphere to be invisible in the line profile and consequently they are also 
not recovered by the inversions. To successfully detect Stokes linear polarization in this quiet snapshot, a noise level below 10~ 3 5 is 
necessary. 
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1. Introduction 

Measuring the magnetic field of the solar chromosphere is an 
outstanding challenge for observers. The reason for this is three- 
fold: there are not many lines available in the optical spectrum 
with a sufficiently high opacity to place the formation in the 
chromosphere, the polarimetric response is limited, and the few 
lines available have a non-local component to their formation 
(|Socas-Navarro & Trujillo Bueno|[T997l |Pietarila et al.||2006; 



|Manso Sainz & Trujillo Bueno|20iO[|Leenaarts et al.|2012| )T 

The situation is much simpler in the photosphere: to model 
most Fei photospheric lines, it is cus tomary to assume loca l 
thermodynamic equilibrium (LTE, see |Rutten~& Kostik 1982), 
which means that the atom population densities can be com- 
puted directly from the local temperature and electron density 
at each point. Unfortunately, most chromospheric lines require 
a non-LTE (NLTE hereafter) treatment of the radiative transfer 
problem. 

Out of the available chromospheric diagnostic lines, the 
Can infrared (IR) triplet (/18498, /18542, /18662) constitutes a 
good compromise of modeling efforts, polarimetric sensitiv- 
ity and observational requirements; partial redistribution effects 
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are negligible ( [Uiten broek 1989 ) and almost the entire calcium 
is singly ionized under typical solar chromospheric conditions 
such that time-dependent ionization effects are negligible and 
the statistical equili brium equations can be used to com pute the 
level populations ( Wedemeyer-Bohm & Carlss on[2011| ). 
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The Ca n IR triplet lines have relatively low effective Lande 
factors (g 8498 = 1.07, g 854 2 = 1.10, g 8662 = 0.87) and broad 
Gaussian line cores compared to photospheric lines. This results 
in low polarization signals compared to those from photospheric 
lines. However, the triplet lines have the advantage that they 
can be observed over the entire solar disk (as opposed to the 
He i /U0830 multiplet, which is visible only in active patches) 
and provide information also on the plasma temperature. 

The development of inversion techniques has improved our 
capabilities to infer physical quantities from the chromosphere 
(see^ujillo Bueno 2010] [Asensio Ramos et aL|2012| and ref- 
erences therein). The first full-Stokes NLTE inversions of the 
Can infrared lines were carried out by |Socas-Navarro et aL 
(2000). Their approach is ideal for analyzing active region and 
network observations where the magnetic field is strong. To an- 
alyze quiet-Sun observations, very high sensitivity is required, 
especially in Stokes Q and U, where signals are expected to be 
very low. 
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This study aims to improve our understanding of chromo- 
spheric magnetic field diagnostics. To this end we will use 3D 
simulations to study NLTE inversions of quiet-Sun observations 
to quantify their reliability. We will also study the detectabil- 
ity of Zeeman-induced polarization; instrumental limitations and 
requirements. The goal is to investigate the polarimetric proper- 
ties of the Ca n IR lines and the regime in which inversions yield 
realistic results. 

A 3D MHD simulation is used to compute synthetic full- 
Stokes profiles with a realistic 3D solution of the NLTE problem. 
The resulting synthetic profiles are inverted assuming plane- 
parallel geometry on a pixel-by -pixel basis. Realistic values of 
noise and instrumental degradation are applied to the simulated 
observations before carrying out the inversions. 

In real observations, Stokes Q and U are produced by 
the combined action of scattering polarization and transverse 
Zeeman effec t when 10 G < B < 100 G ( [Manso Sainz & Trujillo] 
Bue no|20T0| ). Within the mentioned range, Hanle signals are ex- 



pected to have Q and U amplitudes < 10" 4 at solar disk center 
and < 10" 3 close to the limb. In active regions the Zeeman con- 
tribution would be dominant. We emphasize that this study only 
accounts for Zeeman-induced polarization, a simplification that 
is needed given the computational demands of the problem (see 
ICarlin et al.|20l2) . 

The 3D simulation, the radiative transfer calculations and the 
simulated observations are explained in Sec t. [2] the detectability 
of the Stokes signal is discussed in Sect. [3] and the inversion 
code and the inversion results are discussed in Sect. [4] The con- 
clusions are given in Sect. [5] 

2. Simulated observations 

2.1. 3D MHD model 

We used a snapshot from a 3D radiation-MHD simulation of 
the quiet- Sun that extends from the upper convection zone to 
the low corona. The simulation was carried out with the Oslo 
stagger code (OSC, Hansteen et al.||2007| ) which includes non- 
gray radiative losses using multi-group opacities with scatter- 
ing, NLTE radiative losses in the chromosphere, optically thin 
radiative losses in the corona and thermal conduction along 
the magnetic field lines. The simulation domain consists of a 
grid of 256 x 128 x 160 points covering a physical range of 
16.6 x 8.3 x 15.5 Mm, extending from the upper convection 
zone to the lower corona (from 1.5 Mm below to 14 Mm above 
average optical depth unity at 5,000 A). The simulation has an 
average magnetic field strength of 150 G in the photos phere, a 
value consistent with the mean field strength inferred by Trujillo 
Bue no et al.| ( |2004| ) via the Hanle effect. The magnetic field 
is mostly concentrated in intergranular lanes at photospheric 
heights, whereas in the chromosphere it tends to fill up more 
space. At z = 1000 km, the average magnetic field strength is 
79 G. 

Fig. [T] illustrates the vertical component of B in the 
photosphere (bottom) and in the low chromosphere (top). The 
simulation snapsho t us ed in this work is the same as in |Leenaarts 
letaLl ( [20091 and |Rutten et aE1 ( [20TT| ). 



2.2. Radiative transfer 

Our synthetic full Stokes profiles are computed in two 
steps. First, population dens ities are calculated with Multi3D 
(Leenaarts & Carlsson 2009), which allows one to evaluate the 
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Fig. 1. Vertical component of the magnetic field, B z , at horizontal 
slices through the 3D model atmosphere at z=1000 km (top) and 
z= km (bottom). 



three-dimensional radiation field (3D NLTE hereafter). The Can 
atom model used in this work consists of five bound levels plus 
a continuum. We calculate the departure coefficients, defined as 



if: 



(1) 



where ni is the population density of an atomic level i computed 
in NLTE and n* is the equivalent LTE population density. 

In the second step, the same atmospheric model is used 
to compute the LTE solution with the non-LTE inversion code 
based on the Lorien engine (Nicole, [Socas-Navarro et al.|2Q12j ). 
The departure coefficients obtained previously with Multi3D are 
applied to the LTE populations. Then, each column is resampled 
to a grid that is optimized for the radiative transfer computa- 
tions in Can lines: for each column the model is truncated at the 
height where the temperature is higher than 50,000 K. The grid 
points are then distributed according to gradients in temperature, 
density and opacity, placing more points where the gradients are 
steep. The temperature, electron density, mass density, velocity, 
and departure coefficients are interpolated to the new grid using 
an accurate interpolation algorithm based on Hermitian splines 

<Hill|1982). 

Leenaarts et al.| ( |2009| ) pointed out that the Gaussian core 



of the synthetic Ca n /18542 A Stokes / profile is narrower and 
steeper than those from observations, probably from the lack of 
small scale random motions in the model. Assuming Zeeman 
splitting in a weak-field regime, the shape of Stokes Q and V can 
be expressed as a function of the first derivative of the intensity 
( Landi Pegl'Innocenti & Landolfi|2004| ). 



1 



Q(A) = — AAiG— . , 
^ 4 B rf \dA 



dl 



V(A) 



-AAbS cos 6 



(2) 
(3) 



where AA B represents the Zeeman splitting of the line, g is the 
Lande factor, is the inclination, G depends on the quantum 
numbers of the line and 7] is the Voigt profile. According to 
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Wavelength [A] 

Fig. 2. Spatially-averaged profile without micro-turbulence 
(black- solid), with micro- turbulence (gray- solid) and the FTS at- 
las (gray-dashed). 



Eqs.[2]and[3]the amplitude of Stokes Q and V is enhanced if the 
slope of Stokes / is steep. Thus, steeper Stokes / profiles produce 
more strongly peaked Stokes Q, U and V profiles. Therefore, our 
Stokes 2, U and V spectra would be unrealistically strong and 
narrow, which would lead to an overestimation of the effect of 
spectral smearing and an underestimation of the effect of photon 
noise. 

To make more realistic predictions, we introduced an artifi- 
cial micro-turbulence in the model. Without micro-turbulence, 
the full width at half maximum (FWHM) of the chromospheric 
core, measured in our spatially resolved spectra is about 1 80 mA. 
Corresponding measurements from quiet- Sun line profiles were 
carried out by |Cauzzi et al. ( 2009 ), who found values in the range 
450-550 mA. This means that our profiles are of the order of a 
factor 2.5 narrower than observed, assuming that instrumental 
degradation has a negligible effect on the FWHM of the pro- 
files. By artificially introducing micro-turbulence of 3 km s" 1 the 
width of our profiles is increased up to ~ 400 mA. A compar- 
ison between the spatially averaged spectrum with and without 
micro-turbulence is given in Fig. [2] 

After adding the micro-turbulence the formal solution of the 
polarized radiative transfer equation was computed with Nicole 
using a Hermitian method, developed by Bellot Rubio et al. 
( 1998| ). This two-step process from computing the Stokes pro- 
files is necessary because Multi3D cannot compute polarized 
radiative transfer, and Nicole cannot compute 3D radiative trans- 
fer. Only by combining the two codes we obtain full Stokes pro- 
files computed with 3D radiative transfer that serve as our syn- 
thetic observations for the inversion. 

Our synt hetic intensity pr ofiles are near-identical to the re- 
sults of Leenaarts et al. (2009 ), with small remaining differences 
due to small differences in the radiative transfer codes. In Fig. [3] 
monochromatic images at AA = 14.9, 0.84, 0.24, 0.16 and 0.00 
A from the core of the line illustrate the vast formation range of 
the Can /18542 A line, which covers the photosphere and part 
of the chromosphere. The polarization profiles are induced by 
Zeeman splitting. The Stokes Q and U profiles show amplitudes 
about 10" 4 (relative to the continuum intensity), peaking around 
10" 3 at locations where the magnetic field is (relatively) strong 
and mostly horizontally oriented. The Stokes V signal is one or- 
der of magnitude stronger than Stokes Q and U, usually within 
the range [10~ 3 , 10" 2 ], and shows extended areas with the same 
polarity. 

The last column in Fig. [3] illustrates for each wavelength 
the height where the monochromatic optical depth equals one, 
Z(r y =l). In the wings of the line, the height of the r y =l surface 



is dominated by the temperature sensitivity of the opacity. This 
gives rise to granular "hills" with a higher r y =l surface than in 
the cooler intergranular lanes. Additionally, plasma evacuation 
in magnetic elements reduc es the opacity, allo wing to see deeper 
into the atmosphere (e.g. |Keller et al.| 2004} . At A/l=0.837 A 
we sample the inverse granulation layer where the temperature 
contrast is lower and most of the variation in the height of the 
r y =l surface is caused by the lower plasma density where the 
magnetic field is strong. Closer to line center there is a strong 
variation in the r y =l surface caused by strong temperature and 
density variations and Doppler shifts of the line core into or out 
of the passband. 



3. Observability of Stokes signals 

This study is partially motivated by the instrumental re- 
quirements of future telescopes. In this context, Fabry-Perot 
Interferometers (FPI) have become very popular because they 
allow a balanced trade-off between cadence, spatial resolution 
and spectral resolution. We analyze here the combined effects 
of limited spectral resolution and noise in our synthetic observa- 
tions. 

We assumed a Gaussian spectral filter transmission and ad- 
ditive photon noise produced by the finite photon flux. The in- 
strumental point spread function (PSF) is characterized by its 
FWHM, and hereafter we refer to this parameter when instru- 
mental spectral smearing is discussed. The additive noise is char- 
acterized in terms of its standard deviation cr relative to the con- 
tinuum intensity. 

Our results are illustrated in Figs. [4 and [5] at +205 mA from 
the line core, where the red peak of polarization is, on average, 
located. 

In Fig. [4] we show the Stokes Q signal in a subfield of the 
synthetic image for various combinations of spectral degrada- 
tion and noise. The spectral degradation spans the range [0, 100] 
mA, in 25 mA steps, the noise levels are cr = 0, 10" 4 , 10" 3 ' 5 
and 10" 3 . The noise value and spectral degradation are indicated 
inside each panel. 

In the top row of Fig. [4] only spectral degradation is applied. 
Even though the amplitude of Stokes Q is decreased by a factor 
2 between FWHM = mA and FWHM = 100 mA, the main 
features are still visible at maximum spectral degradation. 

The middle rows (panels (f) - (o)) correspond to cr = 10" 4 
and cr = 10" 3 ' 5 . Stokes Q features are still above the noise for all 
assumed values of the FWHM. This result is somewhat expected 
because our Stokes Q profiles usually peak within [10~ 4 , 10" 3 ], 
so noise only affects the very weak polarization signals. 

The bottom row of Fig. [4] corresponds to cr = 10" 3 , where 
noise clearly dominates the Stokes Q images. Some regions 
with strong polarization in Q are discernible in panels (p) 
and (q) (FWHM < 25 mA) but are barely visible when 
FWHM = 50 mA. Above this value, Stokes Q vanishes be- 
low the noise because of spectral degradation. 

Fig. [5] shows Stokes V images for the same subfield as Fig.[4j 
for the same values of spectral degradation. Only the case corre- 
sponding to cr = 10" 3 is shown because the signal is well above 
the noise in all other cases. 

Above, we have described and discussed several combina- 
tions of noise and spectral resolution. However, for a given tele- 
scope aperture, the signal-to-noise ratio can be derived as a func- 
tion of the instrumental efficiency, pixel size and exposure time. 
To estimate the number of photons that can be detected at the 
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Fig. 3. Synthetic full-Stokes images computed from a 12.45 x 8.3 Mm patch of the model. From left to right, the first 4 columns 
show monochromatic Stokes /, Q, U and V images, respectively. The last column shows the height where the monochromatic optical 
depth equals one. From top to bottom, the wavelength decreases from the far red wing to the core of the line. The wavelength relative 
to line core is indicated on each of the Stokes / panels. The Stokes Q, U and V panels have been scaled independently to enhance 
visibility of weak polarization features; the scale range is indicated within each panel in units of the continuum intensity. The 
tick-mark separation is 2 Mm. The panels have been transposed with respect to the panel orientation of Fig. 1. 
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Fig. 4. Stokes 2 image spectrally degraded with Gaussian PSFs and random noise. From left to right, the spectral degradation 
increases as a function of the FWHM of the Gaussian, from AA = (no spectral smearing) to A A = 100 mA. From top to bottom, 
the noise component increases, from cr = to cr = 10" 3 . The images are scaled to ±0.15% relative to the continuum intensity. 
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Fig. 5. Stokes V image spectrally degraded with Gaussian PSFs. From left to right, the spectral degradation increases as a function 
of the FWHM of the Gaussian, from AA = (no spectral smearing) to AA = 100 mA. A Gaussian random noise component of 
cr = 10" 3 has been added. The images are scaled to ±5% relative to the continuum intensity. 



core of the Ca n /18542 A line, w e computed the continuum in- 



tensity at 8542 A using Mum (see 



Carlsson 



atmospheric model by |Fontenla et al. (2006). The intensity ra- 
tio between the continuum and the core of the line (In) is de- 
rived from the Fourier Transform Spectrometer at the McMath- 
Pierce Telescope (hereafter called the FTS atlas) of Bra ult &| 



|Neckel| ( [T987] ), resulting in I = Ii c /I CO nt = 0.18. For a tele- 
1986| and the FALC scope aperture of 1.5 m, a spectral resolution of R = 200000 
(42 mA at /18542), a pixel size of 0.1 " and an overall system 
efficiency of 0.1 (telecope-instrument-detector), one can detect 
n = 2.66 x 10 6 photons per second. This corresponds to a pho- 
ton noise level of 6.1 x 10" 4 . 
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Fig. 6. Stokes Q image spectrally degraded with Gaussian PSFs and random noise. From left to right, the spectral degradation 
increases as a function of the FWHM of the Gaussian, from AA = 25 mA to AA = 100 mA. From top to bottom, the integration time 
is oo, 16 s, 4 s and 1 s, respectively. The images are scaled to ±0.15% relative to the continuum intensity. For a fixed exposure time 
the noise in each image increases with decreasing width of the instrumental profile. 



Fig. [6] illustrates similar results for several combinations of 
exposure times (oo s, 16 s, 4 s and 1 s) and spectral degrada- 
tion (0 mA, 25 mA, 50 mA and 100 mA). One can see the dif- 
ficulties that arise from having a too high spectral resolution: 
as the instrumental transmission profile becomes narrower, po- 
larization features become more prominent but the noise is also 
significantly higher. By increasing the exposure time, noise is 
decreased, but as we explain below, image degradation can ap- 
pear from the rapid chromospheric motions. Additionally, higher 
spectral resolution increases the number of wavelength points 
that are needed to sample the line profile. For clarity, Table [T] 
summarizes the numbers used in Fig. [6] 

The FPI instruments that are currently operating in solar tele- 
scopes can typically achieve a spectral resolution of 45-100 mA 
in the near-infrared, with noise levels of about 10" 3 . Those cases 
correspond to panels (r), (s) and (t) in Fig.|4]and[5] Our results 
suggest that Stokes V detections are possible with current in- 
strumentation. However, the detection of linear polarization in 
the quiet Sun requires sensitivities of at least 10~ 3 5 . 



Table 1. Photon noise levels as a function of the exposure time 
and the spectral resolution. 





25 mA 


50 mA 


75 mA 


100 mA 


16 s 


2.00 x 10" 4 


1.42 x HT 4 


1.16 X 10" 4 


1.00 x 10" 4 


4s 


4.00 x 10" 4 


2.83 x 10" 4 


2.31 x 10" 4 


2.00 x 10" 4 


1 s 


8.02 x 10" 4 


5.67 x 10" 4 


4.63 x 10" 4 


4.00 x 10" 4 



Exposure time in solar observations is limited by the evolu- 
tion time of the Sun. Longer exposures than the evolution time 
scale can produce smearing of the data. Particularly, image re- 
construction assumes that the Sun is static for a given set of short 
exposure images. The vigorously dynamic chromosphere shows 
features moving horizontally within one second ( [van Noort & 
Rou ppe van der Voort|2006| ). Additionally, FPI instruments can 
only observe one wavelength at the time, a limiting factor for 
the wavelength coverage given that all acquisitions within one 
line scan should be as simultaneous as possible. Therefore ob- 
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Fig. 7. Stokes / and V profiles from one of the columns from the 
3D model. The instrumental spectral FWHM is 100 mA. The 
gray circles indicate a coarse grid of points located outside of 
the chromospheric range of the line. The black circles sample 
the core of the line every 50 mA. 

serving the chromosphere involves a difficult trade-off between 
wavelength coverage and exposure time. 

There are ways to ameliorate the sampling. The broad Can 
infrared lines have a wide formation range: the fast-moving fea- 
tures observed in the cores contrast with the slow evolution 
of photospheric granulation present in the wings of the lines. 
Additionally, the amplitude of Stokes Q, U and V rapidly de- 
creases toward the wings of the line. 

A strategy to optimize the polarization measurements and 
achieve proper time sampling would be to place a fine grid of 
points in the chromospheric core and a coarser grid in the wings, 
where slower variations are expected. Because the photosphere 
evolves on a longer time scale than the chromosphere, exposures 
in the wings do not have to be recorded as frequently as in the 
core. Fig. [7] illustrates this scheme, with a fine grid of points 
(black circles) distributed in the core every 50 mA and a coarser 
grid in the wings (gray circles) placed every 500 mA. 

4. Inversions 

We used the code Nicole for the inversions presented in this sec- 
tion. It solves the NLTE problem assuming plane-parallel geom- 
etry, isotropic scattering and complete f requency redistribution, 
using the strategy described by |Socas-Navarro & T rujillo B ueno| 
( |1997| ). All Zeeman sublevels originating from a given atomic 
level are assumed to be equally populated, discarding any quan- 
tum interference between them, as proposed by Trujil lo Bueno| 



& Landi Degrinnocenti| ( |1996| ). 



The velocity-free approximation is used for our calcula- 
tions here, which assumes a static atmosphere when the atomic- 
populations are computed. Under these conditions, only half of 
the profile needs to be computed and fewer quadrature angles are 
needed. Once the NLTE populations are calculated, the velocity 
stratification is taken into account to produce the final emerging 
profiles. 
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Fig. 8. Illustrative examples of the fits (solid line) to the simu- 
lated observations (dashed line) at two different locations in the 
LTE inversion of Fig. [I4l Left panels: A good fit, at pixel coor- 
dinates (25,23) in Fig.|14| Right panels: A comparatively poorer 
fit at pixel coordinates (39,34). Upper row: Stokes / (normalized 
to the average continuum intensity). Lower row: Stokes V. 



The inversions are initialized with a guess model that con- 
tains temperature, electron pressure, line-of-sight (l.o.s.) veloc- 
ity, micro-turbulence and the three components of the mag- 
netic field vector. The inversion engine of Nicole is based on 
a Levenberg-Marquardt algorithm that computes the correc- 
tions to a guess model in order to minimize the differences be- 
tween the observed and the synthetic profiles. To maintain phys- 
ical consistency between temperature, electron pressure and gas 
pressure, the model is set to hydrostatic equilibrium after each 
correction. For a given temperature, the gas pressure and the 
electron pressure are obtained using an equation of state and the 
hydrostatic equilibrium approximation. 

The corrections to the guess model were applied at node 
points that are equidistantly distributed over the depth scale of 
the model. The first inversion was initialized with the VAL-C 
model. After the first inversion another five inversions were per- 
formed starting from the VAL-C model where the atmospheric 
quantities have been given a random offset. This minimized the 
odds of the algorithm settling into a local minimum of the x 2 
hyper surf ace. 

Fig. [8] shows the difference in the fit between a pixel that 
converged to the correct solution and another one for which the 
algorithm settled in a secondary minimum of the x 1 hypersur- 
face. The noise created in the resulting image by imperfect fits 
(we will refer to this as inversion noise to avoid confusion with 
the usual photon noise found in the observations) can be mini- 
mized by implementing more sophisticated schemes with mul- 
tiple initializations to ensure that one reaches the absolute min- 
imum. Genetic algorithms (Charbonneau 1995 ) provide a suit- 
able approach to the problem and have be en successfully im - 
plemented in Stokes inversion codes (e.g., |Lagg et aL||2004 



However, such an algorithm dramatically increases the already 
very large amount of computational work required for the so- 
lution of the NLTE inversion problem. For this reason we have 
refrained from implementing anything beyond the already men- 
tioned five attempts with randomized initializations. 
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Table 2. Summary of node points used during each cycle of the 
inversion. 



Physical parameter 


Nodes in cycle 1 


Nodes in cycle 2 


Temperature 


6 


10 


l.o. s velocity 


2 


5 


B z 


1 


4 


B x 





1 


By 





1 



The spectral resolution of the instrument, the wavelength 
coverage and sampling and the signal-to-noise ratio set con- 
straints on the amount of information that can be retrieved by 
an inversion. For a given set of these quantities there is a maxi- 
mum number of nodes after which adding more nodes does not 
improve the inversion of the data. Adding more nodes can lead 
instead to artifacts. 

To improve convergence, the inversions were computed in 
two cycles (see Ruiz Cobo & del Toro Iniesta 1992| |Pieta rila 



et al. 



2007| |Socas-Navarro||2011| ). In the first cycle, a limited 



set of nodes was used to obtain a rough solution that was then 
refined during the second cycle with more nodes. A summary of 
the number of nodes used on each cycle is shown in Table [2] 

The four Stokes parameters are weighted differently in the 
inversion. The following ratio was used in the present work: 
W(I,Q,U,V) = [1,20,20,8]. 

The inversions were computed with critically sampled pro- 
files. For spectral resolutions of 50 mA and 100 mA, the spec- 
tral sampling at the chromospheric core are 25 mA and 50 mA, 
respectively. The photospheric line wings are covered with a 
coarse grid of points separated by 1 A in all cases. 

4.1. Non-LTE inversions 

We performed NLTE inversions of our simulated Stokes pro- 
files, which were first convolved with a Gaussian function of 50 
FWHM, without noise. 

The results of these inversions are compared with the orig- 
inal atmosphere in Fig. [9] Each quantity that is shown is the 
average over height between log 10 (T5oo)= -4.5 and the depth 
where the transition region starts at each pixel, on average 
log 10 (T5oo)= -5.52. The inversion code usually fails to repro- 
duce the steep temperature gradient of the transition region, 
therefore it has been excluded from the depth average. In the 3D 
snapshot, the log 10 (T5oo)=-4.5 iso-surface corresponds on aver- 
age to z = 737 ± 83 km; the log 10 (T5oo)=-5.52 iso-surface is 
located on average at z = 1340 ± 365 km. We stress that each 
column has its own run of T5000 with height in the atmosphere. 
The derived quantities can therefore not be directly interpreted 
as the average over height in a horizontal slab. 

The first row in Fig. [9] compares temperatures. The inver- 
sion reproduces most of the high-temperature structure, but fails 
to reproduce the low temperatures. This failure is at least partly 
explained by the effect of 3D radiation, which is not taken into 
account in the inversion, because this assumes each pixel is a ID 
plane-parallel atmosphere. Despite the failure to reproduce low 
temperatures, the average relative error in the inferred temper- 
ature (measured as the difference between model and inversion 
divided by the model) is only 6.7%. The Pearson correlation co- 
efficient for the temperature is r = 0.51. If we exclude cold areas 
(T < 3500 K), the correlation factor becomes r = 0.67. 

The second row of Fig. [9] shows the line-of- sight velocity. 
The inversion is able to recover almost the entire structure orig- 



inally present in the simulation, with a Pearson correlation coef- 
ficient of r = 0.80. The differences between the model and the 
inversion are mainly caused by inversion noise. Although in this 
paper we explore the inversion of one spectral line, it is still pos- 
sible to retrieve information about gradients (mostly on the line- 
of-sight velocity and the magnetic field) because these gradients 
are needed to reproduce the line asymmetries. The extreme case 
of a shock propagating through the atmosphere is discussed be- 
low in more detail. 

The third row shows the vertical component of the magnetic 
field. The inversion reproduces the values in the atmosphere very 
well, with a Pearson correlation coefficient of r = 0.97. The 
fourth and fifth row finally show the horizontal components of 
the magnetic field. The quality of the reconstruction is lower 
than for B z , but still very good, with correlation coefficients 
of r = 0.90 for \B X \ and 0.78 for \B y \. We note the presence 
of some artifacts in the inverted \B X \ and \B y \ maps, such as in 
\B x (x = 6.5, y = 4.5)|, where the inverted field is stronger than in 
the 3D snapshot. These artifacts appear because the inversion is 
computed here with one node in \B X \ and one node in \B y \. The 
presence of strong photospheric fields can pollute the inverted 
values. Unfortunately, increasing the number of nodes also in- 
creases the inversion noise. 



Figure 10 shows probability density plots comparing the in- 
verted quantities with the original data. All inverted quantities 
show a tight correlation with the original data, except for the 
temperature. 

The failure to reproduce low temperature areas is further il- 
lustrated in Fig. [TT] The top panel compares the emerging in- 
tensity profile for a pixel with a cold chromosphere (shown in 
the lower panel) computed in 3D and in ID. In ID, the line- 
core intensity is low, as a consequence of the low temperature. 
In 3D, however, the cold gas is illuminated from the sides, which 
increases the source function and hence the emerging intensity. 
The lower panel shows the temperature in the model atmosphere, 
and the temperature from the inversion. The inversion shows a 
much higher temperature (and hence source function) to repro- 
duce the relatively high 3D line-core intensity. 

The presence of shocks in the atmosphere poses a particu- 
lar problem for the inversion. Owing to the limited number of 
nodes, the inversion will often not be able to recover the true at- 
mospheric structure. Instead it will try to find a much smoother 
structure that nevertheless approximates the emerging intensity. 

We illustrate this with an example from the 3D simulation 
in Fig. 
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the shock is located between log 10 (T5oo) = -5 and 
Iog 10 (r5oo1 = -6, as indicated by the increasing temperature and 
a spike in the vertical velocity. Together these lead to an asym- 
metric line core. The profile minimum is blueshifted with respect 
to the rest line-center frequency and the profile shows an inflec- 
tion point at AA = 0.16 A. 

The profile-minimum intensity is well reproduced by the in- 
version, as is the profile-minimum Doppler shift, but the inflec- 
tion point is not recovered. That the profile-minimum intensity 
is recovered means that the inferred temperature follows the true 
temperature in the shock region. The atmospheric velocity is 
not accurately recovered: the limited number of nodes forces a 
smooth inferred velocity field, which yields the same profile- 
minimum Doppler shift as the forward calculation, but has a 
lower maximum velocity. The inferred velocity has a smooth 
transition from 2 to kms" 1 between log 10 (T5oo) = -5 and 
log 10 (T5oo) = -4, whereas the atmosphere has a much steeper 
velocity gradient. The lack of this steep gradient causes the lack 
of the inflection point in the inverted profile. 
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Fig. 9. Full field-of-view comparison between physical quantities in the chromosphere as determined from the NLTE inversion (left) 
and the original 3D snapshot (right). From top to bottom: temperature, line-of-sight velocity, B z , \B X \ and \B y \. The line profiles were 
degraded with an instrumental profile of FWHM = 50 mA before carrying out the inversion. In the velocity panel, upflows have 
positive sign. 



Fig. [13] shows the effect of spectral resolution and photon hand mosaics shows that the quality of the inversion is more 

noise on the quality of the inversions. Because of the large affected by noise levels than by the spectral degradation consid- 

amount of computational work, we restricted the inversions to ered in our inversions. For a given noise level, the results ob- 

an area of 4x2.5 Mm 2 with its lower left corner at (x, y) = (6, 3) tained for a spectral resolution of 50 mA are slightly better than 

Mm in Fig. [9] We have introduced Gaussian noise of cr = at 100 mA. 
10" 4 , 10" 3 ' 5 and 10" 3 . Comparison of the left-hand and right- 
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Fig. 10. Probability density plots of the inverted quantities vs. the 3D simulation quantities. From left to right: temperature, V z , B z , 
\B X \ and \B y \. The Pearson correlation coefficient (r) is indicated in the top-left corner of each panel. 
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Fig. 11. Upper panel: Vertically emerging intensity profile for 
a pixel with a low chromospheric temperature computed with 
Multi3D using 3D radiative transfer (solid-gray), the profile for 
the same pixel, but assuming ID plane-parallel geometry (solid- 
red) and the best fit from an inversion (dashed). Lower panel: 
The temperature structure for the same pixel (solid-gray) and 
the temperature structure inferred from the inversion (dashed). 



The longitudinal magnetic field (B z ) is the most accurately 
retrieved parameter from the inversions. The Stokes V profiles 
usually peak well above the 10" 3 level. The reconstructed longi- 
tudinal field becomes increasingly noisy with increasing photon 
noise, but the overall structure remains the same. 

The maximum synthetic Stokes Q and U signal is typically 
only 5xl0~ 4 . Reconstruction of the horizontal components of 
the magnetic field is therefore much more sensitive to the noise 
level. At 10" 4 noise, the inversion is able to recover most of the 
structure present in the atmosphere, more successfully when the 
FWHM = 50 mA. When noise is increased to cr = 10" 3 5 , most 
of the details are lost, but areas with horizontal field larger than 
~100 G are still discernible. Finally, when the noise is of the 
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Fig. 12. Upper panel: Intensity profile calculated from a col- 
umn of the model where a shock is present (solid line). The 
dashed line corresponds to the inversion result. Middle panel: 
Temperature stratification in the original model (solid) and the 
inversion result (dashed). Bottom: Line-of- sight velocity strat- 
ification in the original model (solid) and the inversion result 
(dashed). The shock is located between the two gray dash-dotted 
lines. 
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order of 10 3 , the inversion is unable to recover the transverse 
magnetic field components. 

4.2. Velocity and magnetic field In LTE 

The intensity absorption pro file of a given spectral line is given 
by the expression (see, e.g. Landi Degl ' Inn ocentij 1 992 ) 



ki(A, T, B, vios) = k c {T) + ^(^ n h n u ) 

<f> b (A, B, Vios) + 0rW, B, Vios) /t om 

+ (1 + cos 20) 



> p (A,B, vi os ) sin 20 
(4) 



where kq is the continuum absorption coefficient, rjo is the line 
strength, <p p , fa and (p r are the n, blue and red profile compo- 
nents, respectively. Eq. [4] explicitly shows the dependencies on 
the temperature, T (and other thermodynamical parameters such 
as density or pressure), magnetic field vector B and the line-of- 
sight velocity vi os . For the continuum absorption coefficient kq 
we have neglected the wavelength dependence because the con- 
tinuum opacity varies with wavelength much more slowly than 
the line opacity. The atomic level populations of the lower and 
upper levels of the transition considered, n\ and n u , obviously de- 
pend on T as well. Even in conditions far from LTE, the atomic 
populations are only weakly dependent on B and vi os , unless the 
velocity gradients are very steep. 

In Eq. [4] there is a convenient separation of the atmospheric 
parameters (the unknowns of our inversions). The dependence 
on the thermodynamical parameters T is contained exclusively 
in kc and rjo, whereas the dependencies on the magnetic field B 
and the velocity vi os are in the term enclosed in square brackets. 
We can then see that the wavelength dependence of the profile 
(which is contained in the terms within the brackets) depends 
only on the values of B and Vi os . The same happens with the 
equivalent absorption profiles for the rest of the Stokes parame- 
ters kq, ku, Ky, as well as the emission and the anomalous dis- 
persion profiles. They can all be separated into a term containing 
the atomic populations and the thermodynamic properties and 
another term containing the wavelength profiles with the depen- 
dence on B and vi os . 

This separation of variables is fortunate because if we are 
not interested in the temperature but only in the magnetic field 
and velocity, it is possible to devise a much simpler inversion for 
these quantities. This simpler inversion would not require the 
solution of an NLTE problem, and the computational demands 
are much lower. 

Let us consider now one such simplistic procedure, namely 
one that fits the observed profiles computing the atomic popu- 
lations in LTE. With the inversion we obtain a run of thermo- 
dynamic parameters T* which is, in general, different from the 
real T in the solar atmosphere. From that T* we compute LTE 
populations rf { and n* u . 

When the observations are fitted by the procedure, the source 
function and the opacity are adequately reproduced. Therefore, 
the populations obtained are approximately correct (rtf n\ and 
n* u ^ n u ), even though the temperatures needed to yield those 
values are drastically different in LTE and NLTE. The wave- 
length dependence of the observed spectral line depends entirely 
on B and v/ 0J , whose effects on the profile are the same regard- 
less of whether the populations are computed in LTE or NLTE. 
In summary, with our simplistic LTE fitting procedure we would 
recover B and v/^with the same degree of accuracy as with an 
NLTE inversion, and only T would be wrong. 



With the same reasoning, we argue that questioning the va- 
lidity of the assumptions that we have made on the radiative 
transfer (such as hydrostatic equilibrium, 1.5D atmosphere or 
the assumption of statistical equilibrium) would only be relevant 
in determining T but not B or vi os . 

This approach is actually similar to what is often performed 
in Milne-Eddington inversions for the photosphere, in which 
only B and vi os are interpreted as physical parameters, even 
though several other parameters are involved in the fit. 

Figure [T4| depicts the result of this LTE inversion of the syn- 
thetic Stokes profiles to determine B and vi os in the same re- 
gion of the simulation as for Fig. 13 The inversions in Fig. 14 
correspond to cr = 10" 4 and a spectral resolution of 50 mA. 
Comparison of the inferred B z and vi os shows that the inversion 
performs reasonably well in reproducing these quantities. The 
noise comes primarily from inversions that failed to converge. 
The Pearson correlation coefficients are r = 0.93 for B z and 
r = 0.84 for v hs . 

For B x and B y the results appear more noisy than for B z , but 
the inversion more or less correctly reproduces the patches of 
strong field on the right side of the panel. The inversion erro- 
neously infers patches of horizontal field where there are none 
in the atmosphere, such as at (x,y) = (12, 28). The Pearson cor- 
relation coefficients are r = 0.90 and r = 0.69 for B x and B y , 
respectively. The lower quality of the inversion compared to the 
B z case is due to the Stokes Q and U signals being weaker than 
Stokes V. 

In a realistic case, the data would be dominated by the ob- 
servational noise, except in specific cases such as strong fields in 
the vicinity of sunspots. 



5. Discussion and conclusions 

We have computed full Stokes synthetic line profiles using a 3D 
radiation-MHD simulation of the outer atmosphere of the Sun. 
A full 3D solution of the NLTE radiative transfer problem was 
performed using Multi3D. The width of the core of the synthetic 
Can intensity profiles is approximately a factor 2.5 lower than 
observed. This would lead to artificially strong Stokes Q, U and 
V signal. To remedy this, external micro-turbulence was added, 
which brought the synthetic line widths closer to the observed 
ones. 

The resulting line profiles were degraded by various amounts 
of spectral smearing with various levels of noise added. They 
were then inverted using the inversion code Nicole, which solves 
the NLTE radiative transfer problem assuming plane-parallel ge- 
ometry for each pixel (1.5D). The inversions were performed as- 
suming that the polarization is produced by the Zeeman effect 
only. This assumption is reasonable for solar regions with mod- 
erate to strong activity, which is the most likely target for an 
application of Nicole. 

We find here that the NLTE inversion can reliably recover the 
original line-of- sight magnetic field and velocity, with the noise 
levels and spectral resolutions considered. Temperatures are less 
reliably inverted. High-temperature areas are well reproduced, 
but the inversion is less sensitive to cold areas because it cannot 
properly treat the 3D effect caused by incoming radiation from 
hotter neighboring pixels. 

The quality of the NLTE inversion for the horizontal field 
strongly depends on the amount of noise present in the Stokes 
profiles, and less so on the spectral resolution. Only with a noise 
level below cr = 10" 3 5 , is the inversion able to recover the fea- 
tures present in the model atmosphere. 
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We also tested the reliability of a simpler and faster LTE 
inversion, which can recover the line-of-sight velocity and the 
magnetic field, but not the temperature. We showed that this 
inversion recovers B z and vi os quite well, but yields less reli- 
able results for the horizontal components of the magnetic field. 
Nevertheless, LTE inversion is a viable and fast strategy if one is 
interested in the vertical velocity and magnetic field only. 

Our results strongly suggest that to be able to detect and 
invert Stokes Q and U signals in the quiet- Sun, observations 
should have a noise level better than cr = 10" 3,5 . Instrumental 
filter widths below 50 mA do not improve the quality of the in- 
version. 

Our work leads to a better understanding of the information 
contained in chromospheric Ca n observations and what can be 
learned from their inversions. We showed that a NLTE inversion 
can reliably recover the vertical magnetic field and velocity at 
noise levels of cr = 10" 3 from a radiation-MHD simulation of 
quiet Sun, and that the horizontal components can be recovered 
if the noise level is below cr = 10" 3,5 . Our results also strongly 
suggest that if the inversion techniques are applied to observa- 
tions of more active regions on the Sun, which produce stronger 
polarization signals, the inferred quantities will be fairly accu- 
rate. 
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Fig. 13. Illustration of the effect of spectral resolution and photon noise on the quality of the inversions. The left-hand mosaic shows 
the results of the inversion after spectral degradation of the line profiles with a Gaussian instrumental profile with FWHM = 50 mA, 
the right-hand panel after smearing with a profile with a FWHM of 100 mA. Within each mosaic, from left to right, the panels show 
B z (x, v), \B x (x,y)\, \B y (x,y)\. Rows, from top to bottom, show the magnetic field from the 3D snapshot, followed by the inversion 
results assuming random photon noise of cr - 0, cr = 10" 4 , cr = 10" 3 5 and cr - 10" 3 . 
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Fig. 14. Results from an inversion carried out in LTE. The top row illustrates quantities from the 3D snapshot, the middle row the 
results of the NLTE inversion and the bottom row the results of the LTE inversion. In the velocity panel, upflows have positive sign. 
The inversions were computed assuming a spectral resolution of 50 mA and random noise with cr = 10" 4 . 



14 



